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MANAGEMENT  UITTREKSEL 


Titel  :  IMPACT  ON  A  STEEL  SONAR  DOME 

Auteurs  :  dr.ir.  R.H.B.  Fey,  ing.  JJ.  Uwland 

Datum  :  20  maart  1998 

Opdrachtnr.  :  A97/KM/154 

Rapportnr.  :  98-CMC-R0328 

Aanleiding  voor  het  onderzoek 

In  het  verleden  zijn  de  sonar  domes  van  de  S-fregatten  uitgevoerd  in  staal.  Over  een  periode 
van  meer  dan  15  jaren  zijn  er  schadegevallen  geweest  aan  de  stalen  domes,  maar  die  hebben 
nooit  geleid  tot  disfunctioneren  van  de  sonar,  Vanwege  betere  transmissie-eigenschappen  zijn 
de  domes  van  de  M-fregatten  gemaakt  van  composiet.  Tot  op  heden  is  een  van  de  composieten 
sonar  domes  totaal  verloren  gegaan,  terwijl  er  ook  een  door  impact  zware  schade  heeft  opgelo- 
pen. 

Doel  van  het  onderzoek 

Voor  een  representatief  schadegeval  aan  een  stalen  sonar  dome  zal  worden  geanalyseerd  wat  de 
contactkracht  tussen  het  object  en  de  dome  is  geweest  tijdens  de  impact.  Bovendien  zal  worden 
berekend  welke  hoeveelheid  energie  door  de  stalen  dome  is  geabsorbeerd.  Met  deze  gegevens 
zal  in  een  later  stadium  getracht  worden  specificaties  te  formuleren  voor  een  verbeterde  com¬ 
posieten  dome. 

Korte  omschrijving  van  het  onderzoek 

Ten  einde  de  contactkracht  en  de  geabsorbeerde  energie  tijdens  de  impact  te  berekenen  is  een 
eindige-elementenanalyse  uitgevoerd.  De  moeilijkheid  bij  de  analyse  was  met  name  het  feit  dat 
gegevens  over  de  impactor,  die  de  schade  aan  de  stalen  dome  veroorzaakt  had,  ontbraken. 
Dientengevolge  zijn  aannames  over  vorm,  grootte  en  impactrichting  van  de  impactor  noodza- 
kelijk  geweest. 

Conclusies 

Het  berekende  deformatiepatroon  komt  goed  overeen  met  het  deformatiepatroon  zichtbaar  op 
de  beschikbare  foto’s.  Bij  een  penetratiediepte  van  0.2  m  is  de  berekende  contactkracht  230  kN 
en  de  bijbehorende  geabsorbeerde  deformatie-energie  is  32.5  kJ. 
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1.  Introduction 


In  the  1980’s  sonar  domes  were  constructed  out  of  steel.  Nowadays  they  are  made  of  compos¬ 
ite  material  because  composite  material  disturbs  the  sonar  signal  less  than  steel.  Recently,  one 
composite  dome  was  heavily  damaged  and  one  composite  dome  was  even  completely  lost  be¬ 
cause  of  impact  with  unknown  objects.  The  steel  domes  were  also  regularly  damaged  but  the 
damage  was  always  relatively  local  and  the  sonar  dome  did  never  fail. 

The  goal  of  this  research  is  primarily  to  derive  a  criterion  for  impact  resistance  which  will  be 
used  in  the  design  of  new  composite  domes.  In  order  to  achieve  this,  we  will  study  a  case  of 
impact  on  steel  dome  no.  3,  which  was  damaged  by  an  unknown  object,  see  Figure  1-Figure  8 
(photos  dd.  July  26,  1984).  This  damage  is  thought  to  be  representative  for  impacts  under  oper¬ 
ating  conditions.  More  specifically,  the  amount  of  deformation  energy  which  is  absorbed  dur¬ 
ing  the  impact  by  the  steel  dome  will  be  estimated. 

Because  the  object,  which  hit  the  dome,  is  unknown,  it  is  necessary  to  make  assumptions  about 
the  size  and  the  shape  of  the  object,  the  stiffness  and  mass  of  the  object,  its  speed  during  im¬ 
pact  and  the  direction  of  the  impact. 

The  calculations  presented  in  this  report  were  carried  out  with  the  explicit  finite  element  pack¬ 
age  LS-DYNA3D  [1,2]. 
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2.  Description  of  the  model 


2.1  Geometry 


Photos  of  the  undamaged  dome  are  given  in  Figure  1 -Figure  3.  The  thickness  of  the  hull  of  the 
dome  is  4  mm.  There  are  primary  frames  (160  mm  x  3  mm)  and  secondary  frames  (60  mm  x  3 
mm)  on  the  inside  of  the  dome  in  horizontal  and  vertical  direction.  On  top  of  the  primary 
frames  are  tubes  with  a  radius  of  7.5  mm.  The  horizontal  frame  distance  between  the  vertical 
frames  is  maximal  375  mm.  The  vertical  frame  distance  between  the  horizontal  frames  is  120 
mm. 


The  geometry  of  the  mesh  of  the  dome,  which  will  be  described  in  more  detail  in  the  next  sec¬ 
tion,  is  based  on  drawing  Z  01.10  (August  1976)  entitled  ‘Spantenlijst  (vorm  van  de  boeg- 
dome)’,  drawing  G  44.6  (January  1975)  entitled  ‘Ensemble  du  bulbe’  and  the  telefax  dd.  Sep¬ 
tember  9“^  1997  from  Rijkswerf  Den  Helder  to  MARTECH-MMT  and  the  visit  of  TNO  Centre 
for  Mechanical  Engineering  to  Rijkswerf  Den  Helder  dd.  November  19‘*'  1997. 

As  stated  before,  the  object,  which  hit  the  dome,  is  unknown.  In  order  to  get  an  impression  of 
the  geometry  of  the  object  which  hit  the  dome,  the  photos  of  the  damage  to  dome  no.  3  have 
been  examined  (Figure  4-Figure  8).  On  the  outside  the  damaged  area  has  approximately  a  cir¬ 
cular  shape  with  a  radius  of  0.4  m.  In  this  circle  a  right-angled  dent  is  visible,  see  Figure  4  and 
Figure  5.  Several  hypotheses  can  be  given  for  explaining  the  right-angled  dent: 

1.  On  the  inside  of  the  dome  two  large  stiffeners  are  visible  in  the  damaged  area:  a  longitudi¬ 
nal  stiffener  and  a  transversal  stiffener,  see  Figure  6. 

2.  The  direction  of  impact  is  of  influence. 

3.  The  object  which  hit  the  dome  had  sharp  edges. 

4.  A  combination  of  the  three  hypotheses  given  above. 

On  the  inside  of  the  dome  (Figure  6)  we  see  tilting  of  the  top  horizontal  primary  stiffener  over 
three  vertical  secondary  stiffener  distances.  Also  local  buckling  of  secondary  stiffeners  is  visi¬ 
ble. 

We  will  assume  that  the  object  which  hits  the  dome  is  a  sphere  (because  of  the  circular  shape 
of  the  damaged  area)  with  a  radius  of  0.2  m.  The  radius  of  the  sphere  is  chosen  smaller  than  the 
radius  of  the  damaged  area  (0.4  m)  because  it  is  expected  that  a  plastic  hinge  with  a  circular 
shape  will  move  outward. 
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2.2  Meshing 


Because  the  damage  to  the  dome  is  local,  located  on  one  side  and  relatively  high,  only  the  port 
side  is  modelled.  Actually,  the  dome  visible  on  the  photos  was  hit  on  starboard  side.  This  is  not 
relevant  for  the  analysis,  however,  because  starboard  and  port  side  are  symmetric.  Figure  9 
shows  the  place  on  the  dome,  which  was  damaged  (marked  by  dotted  curve). 

Figure  10-Figure  14  show  the  element  mesh  from  different  view  points.  The  model  consists  of 
11186  shell  elements  and  321  beam  elements  and  has  10735  nodes.  Near  the  location  of  impact 
the  mesh  has  been  refined.  In  this  region  3  shell  elements  are  used  in  the  direction  of  their 
height  for  the  secondary  frames  and  7  shell  elements  for  the  primary  frames.  The  tube  on  top  of 
the  primary  frames  is  modelled  with  beam  elements.  The  minimal  element  length  is  20  mm. 

The  positive  x>axis  points  from  the  aft  to  the  front  of  the  ship.  The  positive  y-axis  points  from 
port  side  to  starboard  side.  The  positive  z-axis  thus  points  from  the  top  to  the  bottom  of  the 
ship. 


2.3  Material  Properties 


The  sonar  dome  is  manufactured  of  AISI  316  L  steel  (Stoff-Nr.  1.4404),  see  Table  1.  The  mate¬ 
rial  model  used  in  the  calculations  is  an  elasto-plastic  model  with  linear  hardening  for  the  shell 
elements  and  an  elasto-plastic  model  without  hardening  for  the  beam  elements. 


Young’s  modulus  [N/m^] 

2.0  10" 

Poisson  ratio  [-] 

0.3 

mass  density  [kg/m^] 

7980 

0.2  %  yield  stress  [N/m^] 

190.0  10‘ 

tangential  stiffness  [N/m^]* 

754.0  10‘ 

Table  1  Material  properties 


*The  tangential  stiffness  is  calculated  here  as  the  difference  between  the  tensile  strength  (490.0 
10®  N/m^)  and  the  0.2%  yield  stress  divided  by  the  failure  strain  (0.4).  The  actual  tangential 
stiffness  could  be  higher  because  it  is  very  well  possible  that  the  maximum  stress  occurs  at  a 
lower  strain,  say  0.2.  Unfortunately,  the  stress-strain  curve  is  not  available  to  verify  this.  This 
implies  that  we  use  a  conservative  approach  here. 

The  impacting  sphere  is  assumed  to  be  rigid. 
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2.4  Boundary  Conditions  and  Initial  Conditions 


The  top  the  dome  is  mounted  to  the  ship  with  bolts.  Here,  displacements  and  rotations  of  the 
dome  are  suppressed.  Because  the  impact  takes  place  near  frame  170  (of  the  ship)  it  is  assumed 
to  be  sufficient  to  model  the  dome  as  far  as  frame  166  (see  drawing  Z  01.10,  August  1976,  en¬ 
titled  ‘Spantenlijst’).  Behind  frame  166  the  dome  is  assumed  to  be  rigid  (the  back  of  the  dome 
lies  against  a  wooden  frame).  The  dome  has  a  very  stiff  double  bottom  which  is  assumed  to  be 
rigid  too. 

On  the  vertical  plane  which  cuts  the  dome  in  half  in  longitudinal  direction,  the  displacements 
in  y-direction  and  the  rotations  around  the  x-axis  and  the  z-axis  are  suppressed. 

The  impacting  rigid  sphere  moves  in  the  horizontal  xy-plane  under  an  angle  of  20  degrees  with 
the  projection  in  the  horizontal  xy-plane  of  the  normal  of  the  dome  surface  at  the  location  of 
impact.  This  projected  normal  makes  an  angle  of  50  degrees  with  the  positive  x-axis.  The  rigid 
sphere  thus  moves  in  the  horizontal  xy-plane  under  an  angle  of  70  degrees  with  the  positive  x- 
axis. 

The  sphere  has  a  constant  velocity  of  10  m/s.  In  fact  the  velocity  of  the  sphere  is  not  of  much 
importance  in  this  analysis,  because  we  aim  for  a  quasi-static  analysis.  This  means  that  the  ve¬ 
locity  of  the  sphere  should  be  small  enough  to  keep  the  amount  of  kinetic  energy  in  the  dome 
during  impact  low.  At  the  same  time  the  velocity  should  be  as  high  as  possible  to  ensure  that 
the  amount  of  cpu-time  needed  for  the  analysis  remains  acceptable. 


2.5  Contact  Algorithm 


The  contact  algorithm  which  is  used  in  the  calculations  is  a  penalty  stiffness  formulation  as  in 
surface  to  surface  contact,  see  [2].  In  the  contact  model  friction  was  not  taken  into  account. 
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3.  Results 


Based  on  Figure  4,  the  real  permanent  penetration  depth  after  unloading  is  estimated  to  be  0.2 
m.  Some  of  the  calculated  results  are  presented  in  time  history  plots.  The  time  can  easily  be 
converted  into  the  penetration  depth  by  multiplying  the  time  with  the  penetration  speed  of  10 
m/s.  This  means  that  in  the  calculations  carried  out  below,  the  penetration  depth  increases  line¬ 
arly  with  respect  to  the  time  and  there  will  be  a  time  dependent  contact  force  between  the  im- 
pactor  and  the  dome:  unloading  is  not  carried  out.  This  means  that  the  permanent  penetration 
depth  (plastic  deformation)  will  be  somewhat  lower  than  the  calculated  penetration  depth 
(elastic  and  plastic  deformation). 


3.1  Deformations 


Figure  15  and  Figure  16  show  the  deformation  of  the  inside  of  the  dome  at  times  t  =  0.02  s  and 
t  =  0.03  s  respectively.  The  deformation  patterns  at  these  two  times  are  not  significantly  dif^ 
ferent.  The  upper  horizontal  primary  stiffener  is  tilted  over  three  vertical  secondary  stiffener 
distances.  Secondary  stiffeners  buckle  locally  (especially  the  horizontal  ones  near  the  fourth 
vertical  secondary  stiffener).  This  deformation  pattern  corresponds  remarkably  well  with  the 
deformation  pattern  visible  in  Figure  6.  The  vertical  primary  stiffener  is  buckled  near  the  lower 
horizontal  primary  stiffener. 

Figure  17-Figure  20  show  the  contour  plots  of  the  displacement  in  x,  y  and  z-direction  and  the 
total  displacement  for  t  =  0.02  s  respectively.  The  damaged  area  is  approximately  circular  with 
a  radius  of  0.4  m,  which  is  the  same  radius  as  calculated  from  Figure  4.  Therefore  the  situation 
at  /  =  0.02  s  will  be  investigated  in  detail.  The  centre  of  the  damaged  area  has  a  maximum  dis¬ 
placement  of  0.2  m. 

Figure  21  shows  the  effective  plastic  strain  at  t  =  0.02  s.  The  maximum  effective  plastic  strain 
is  0.17  and  occurs  in  the  horizontal  secondary  stiffener  just  below  the  region  of  impact. 

In  Figure  7  and  Figure  8  cracks  are  visible  which  are  not  predicted  by  the  model.  The  real 
structure  is  chain  welded.  In  the  model  continuous  welding  is  assumed. 

Figure  22  gives  the  total  displacement  for  t  -  0.03  s. 


3.2  Stresses 


Figure  23-Figure  25  show  the  stresses  in  x,  y  and  z-direction  respectively  for  t  —  0.02  s.  In  the 
region  where  there  is  direct  contact  between  the  sphere  and  the  dome  we  see  tensile  stresses  up 
to  314  MPa.  Just  outside  the  region  of  contact  we  see  compressive  stresses  up  to  -242  MPa. 
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Figure  26  shows  the  Von  Mises  stress  for  t  =  0.02  s.  The  maximum  Von  Mises  stress  is  321 
MPa  and  is  found  in  the  secondary  stiffeners  in  the  region  of  contact. 

3.3  Forces 

In  order  to  check  the  order  of  magnitude  of  the  calculated  forces  two  analytical  formulas  were 
applied:  Roark’s  formula  and  Lupker’s  formula. 

3.3.1  Roark’s  formula 

In  [4]  an  analytical  expression  is  given  for  a  concentrated  transversal  force  acting  on  the 

centre  of  a  flat  circular  elastic  plate  (radius  r )  with  a  simply  supported  edge  as  function  of  the 
transversal  displacement  y  of  the  centre: 

16;cDy  l  +  v 

""  3  +  V 

where: 

Et^ 

12(l-v") 

and  with  V  =  0.3,  E  —  2.0  10*^  N/m^  t  =  0.006  m  (equivalent  thickness  due  to  correction  for 
secondary  stiffeners)  and  r  =  0.4  m  (radius  of  damaged  area). 

For  a  deflection  of  y  =  0.2  m  the  transversal  force  becomes  100  kN. 

3.3.2  Lupker’s  formula 

In  [5]  a  similar  force-deflection  relationship  is  derived  where  a  localised  zone  of  plastic  defor¬ 
mation  and  a  spherical  profile  is  taken  into  account: 

F,  =  I  nM,y{(l-6P,)ky/t}  +  P,ky  / 1] 

where: 

M^  =  a/  lA 

and: 

Pa  =  prl(2a^t) 
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Using  the  values  given  in  the  previous  section  and  (Jq  =  190  10^  Win  (maximum  stress  elas¬ 
tic  perfectly  plastic  material  model),  k  -3  (spherical  profile),  /?  =  0  Wrn  (internal  overpres¬ 
sure),  a  deflection  of  y  =  0.2  m  requires  a  force  =  140  kN. 

For  an  internal  overpressure  of  1  atmosphere  (10^  N/m^)  increases  to  200  kN. 

3,3.3  LS-DYNA3D  calculation 

Figure  27  and  Figure  28  show  the  force  magnitude  and  the  force  in  global  y-direction  as  a 
function  of  time.  From  these  figures  it  becomes  clear  that  the  force  in  global  y-direction  is 
dominant  for  the  force  magnitude.  At  /  =  0.02  s  the  penetration  depth  is  about  0.2  m  and  the 
contact  force  is  230  kN,  which  is  somewhat  higher  than  but  in  the  same  order  of  magnitude  as 
the  forces  following  from  Roark’s  formula  and  Lupker’s  formula.  Probably  it  is  somewhat 
higher  because  the  bending  stiffness  of  the  stiffeners  is  underestimated  in  the  formulas;  moreo¬ 
ver  the  primary  stiffeners  have  been  replaced  by  secondary  stiffeners  in  the  calculation  of  the 
equivalent  thickness. 

The  force  does  not  seem  to  increase  much  more  at  /  =  0.03  s. 


3.4  Energies 


Figure  29-Figure  33  respectively  show  the  external  work,  the  internal  energy,  the  kinetic  en¬ 
ergy,  the  hourglass  energy  and  the  total  energy  of  the  dome  as  a  function  of  time.  The  total  en¬ 
ergy  is  equal  to  the  sum  of  the  internal  energy,  the  kinetic  energy  and  the  hourglass  energy. 
Note  that  the  external  work  should  be  equal  to  and  in  fact  equals  the  total  energy.  Note  further 
that  at  t=0.03s  the  kinetic  energy  is  about  5%  of  the  total  energy  and  the  hourglass  energy  is 
about  2%  of  the  total  energy.  This  means  that  93%  of  the  external  work  is  transformed  in  de¬ 
formation  (internal)  energy  (=summation  of  elastic  and  plastic  energy).  The  deformation  en¬ 
ergy  is  32.5  kJ  for  t  —  0.02  s  and  51  kJ  for  f  =  0.03  s.  At  this  point  we  remind  the  reader  that 
the  approach  we  used  to  calculate  the  tangential  stiffness  is  conservative:  the  actual  tangential 
stiffness  could  be  considerably  higher  (see  page  6).  In  the  latter  case  the  deformation  energy 
will  also  be  higher,  maybe  even  up  to  25%. 

Ai  t  =  0.03  s  the  internal  energy  seems  to  increase  almost  linearly. 

Figure  34  shows  a  contour  plot  of  the  internal  energy  density,  which  makes  clear  to  which  ex¬ 
tent  the  several  parts  in  the  structure  contribute  to  the  internal  energy. 
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4.  Conclusions  and  Recommendations 


The  main  conclusions  are: 

L  The  calculated  deformation  pattern  corresponds  very  well  with  the  deformation  pattern 
found  in  reality.  Both  show  global  tilting  of  the  top  horizontal  primary  frame  over  three 
vertical  secondary  frame  distances.  Also  the  places  at  which  local  buckling  of  secondary 
stiffeners  occurs  correspond  well.  The  radius  of  the  damaged  area  is  in  reality  equal  to  0.4 
m.  In  the  calculation  we  have  this  situation  for  t  =  0.02  s  (penetration  depth  0.2  m).  The 
only  (minor)  difference  between  the  calculation  and  reality  is  that  in  the  last  situation  a 
right-angled  dent  is  visible,  which  is  absent  in  the  calculation. 

2.  Based  on  the  reported  assumptions,  we  find  for  a  penetration  depth  of  0.2  m  a  contact  force 
of  230  kN  and  a  corresponding  deformation  (internal)  energy  of  32.5  kJ.  This  is  20%  of  the 
energy  (165  kJ)  of  an  anchor  (2100  kg)  falling  from  a  height  of  8  m  [3]. 

Recommendations : 

1.  It  would  be  interesting  to  carry  out  a  calculation  for  a  composite  dome  under  the  same  im¬ 
pact  conditions  in  order  to  compare  the  resultant  deformations,  stresses,  forces  and  energies. 
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Figures 


Figure  1  Front  view  of  the  undamaged  dome 
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Figure  2  View  of  starboard  side  of  the  undamaged  dome  (squares  are  numbered  A  to  Dfrom  top  to  bot¬ 
tom  and  I  to  XI  from  left  to  right) 
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Figure  3  Inside  front  view  of  the  undamaged  dome 
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Figure  4  Damaged  surface  has  approximately  a  circular  shape  with  a  right-angled  dent  in  it. 


Figure  5  Outside  view  of  damage 
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Figure  6  Inside  view  of  damage:  tilting  of  horizontal  primary  stiffener  and  local  buckling  of  horizontal 


secondary  stiffeners 


Figure  7  Crack  in  connection  between  hull  and  longitudinal  stiffener 
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Figure  8  Crack  in  stiffener  connection 
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Figure  10  Mesh  of  the  dome  seen  from  the  bottom  in  xy-plane  ( impacting  sphere  in  blue) 
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Figure  13  Mesh  of  the  dome  seen  from  the  inside  near  the  region  of  impact 
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Figure  14  Further  detail  of  the  mesh  of  the  dome  near  the  region  of  impact 
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Figure  15  Deformed  geometry  seen  from  the  inside  at  t-0.02  s 
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Figure  16  Deformed  geometry  seen  from  the  inside  at  t=0,03  s 
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Figure  17 x-displacement  [m]  at  t=0.02  s 
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Figure  18  y-displacement  [m]  at  t-0.02  s 
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Figure  19  z-displacement  [m]  at  t=0.02  s 
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Figure  20  Total  displacement  [m ]  at  t=0.02  s 
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Figure  27  Effective  plastic  strain  [-]  at  t-0.02  s 
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Figure  22  Total  displacement  [m]  at  t=^0.03  s 
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Figure  23  x-stress  [Pa]  at  t=0.02  s 
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Figure  24  y-stress  [Pa]  at  t=0.02  s 
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Figure  25  z-stress  [Pa]  at  t=0.02  s 
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Figure  26  Von  Mises  stress  [Pa]  at  t-0.02  s 
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Figure  27  Total  contact  force  [N]  versus  time  [sj 
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Figure  28  y -force  [N]  versus  time  [s] 
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Figure  29  External  work  [J]  versus  time  [s] 
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Figure  30  Internal  energy  [J]  versus  time  [s] 
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Figure  31  Kinetic  energy  [J]  versus  time  [s] 
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Figure  32  Hourglass  energy  [J]  versus  time  [s] 


2.20E.02 


Total  Energy 


TNO-report 


98-CMC-R0328 


20  March  1998 


SONAR  DOME 

time=  2.00986E-02 
internal  energy  density 

min=  O.OOOE+00  in  element  11617 
max=  4.289E+07  in  element  6885 


fringe  levels 
O.OOOE+OO  > 


8.577E+06> 


1.287E+07  > 


1.715E+07> 


B 


4;289E+07  > 


LS-TAHmiS  936.01  Kch96  I 


Figure  34  Internal  energy  density 
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